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We present an analysis of the quantum state resulting from the dissociation of diatomic molecules 
prepared in a condensate vortex state. The many-body state preserves the rotational symmetry 
of the system in quantum correlated states by having two equally populated components with 
angular momentum adding to unity. A simple two-mode analysis and a full quantum field analysis 
is presented for the case of non-interacting atoms and weak depletion of the molecular condensate. 
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I. INTRODUCTION 

Since the 1995 experiments with the first production of 
atomic Bose-Einstein condensates, degenerate quantum 
gasses have constituted a very active field of research. A 
wide variety of means exists for detection and control of 
the properties of these systems, and the early works have 
been followed by progress on degenerate fermionic sys- 
tems, on mixtures of different species, and on conversion 
between atomic and molecular quantum gasses. 

The coherence properties of the degenerate quantum 
states have been verified implicitly by measurements of 
the response properties of the systems and explicitly by 
the observation of robust interferences £Q and topologi- 
cal structures [2] . Recently, a vortex lattice in a system 
of fermionic atoms was shown to extend over the BCS- 
BEC crossover towards formation of bosonic diatomic 
molecules on the molecular side of a field controlled Fes- 
chbach resonance [3]. 

Many properties of condensates are well described by 
mean field theories and the Gross-Pitaevskii equation, 
but in some cases mean field theories may completely 
fail and totally forbid processes, which occur peacefully 
according to a full quantum state analysis. Important 
examples of such processes are the period-doubling ob- 
served in a shaken lattice [I] and the break-up of a 
moving condensate in a lattice due to four wave mix- 
ing [5] . Both processes are driven by collisions of pairs of 
condensate atoms which emerge in two new momentum 
states. Within a mean field theory, these processes, like 
the equivalent parametric amplification process in opti- 
cal down-conversion and four wave mixing will only get 
initiated if a non- vanishing mean field with the final state 
character is seeded to the solution of the Gross-Pitaevskii 
equation |6J. From a formal perspective, the fact that 
only the sum of the phases of the two final state com- 
ponents and not the individual phases are locked to the 
initial state wave functions causes the fields to remain in 
the vacuum state in order not to break the phase sym- 
metry. If the physical problem is simple, e.g., in the few- 
mode quantum optical problems, one may alternatively 
have recourse to a full quantum many-body theoretical 



analysis. In the context of quantum gasses, we have pre- 
viously [7] studied the "degenerate down-conversion" of a 
molecular condensate by dissociation to a single trapped 
atomic condensate in such a full quantum treatment and 
shown that despite the absence of an atomic condensate 
phase, the atomic component does indeed accumulate in 
a single preferred quantum state. Dissociation into two 
atomic beams have smilariy been proposed as a means for 
generation of Einstein- Podolsky- Rosen correlations [5] . 

Note that in four-wave mixing processes a mean-field 
solution not only breaks the phase symmetry of the field, 
it also breaks the translational invariance of the problem 
as dramatically observed as period doubling in [4], and 
likewise the different populated momentum components 
in [5 will have spatial interference patterns with periods 
exceeding the one of the lattice potential used in the ex- 
periments. In this paper, we consider a process where 
the symmetry breaking is similarly spectacular, namely 
the down conversion by photo dissociation of a molecu- 
lar condensates which is prepared in a single vortex state. 
In a cylindrically shaped trap, a vortex state is a topo- 
logically stable state of a quantum gas with a vanishing 
density along the cylinder axis and a phase which changes 
by 2ir as one follows a closed loop around the axis. The 
rotational symmetry around the z axis is a fundamen- 
tal property of the many-body Hamiltonian and both 
mean field solutions and more elaborate theories must 
respect this symmetry. Microscopically, the vortex solu- 
tion is consistent with a Hartree product state descrip- 
tion with a product of single particle M — 1 eigenstates 
of the azimuthal angular momentum. In this state, ev- 
ery molecule has a unit angular momentum around the 
condensate axis, which cannot be transferred to a pair 
of atoms populating the same (m = 1/2 ?) quantum 
state. Conversely, any mean field solution for the atomic 
system will necessarily correspond to an integer angular 
momentum per atom, and would hence imply an even 
angular momentum of every pair, i.e, of the molecule. 

The molecules must dissociate to at least two different 
atomic states with azimuthal quantum numbers which 
add to the molecular value, m + m! = M. The easiest 
situation to deal with is one where only two states with 
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angular momentum m — and m! = 1 get populated. 
We shall address this case in Sec. II. In particular we shall 
discuss to which extent a strongly number correlated 
quantum state of this kind state is distinguishable from 
a state which macroscopically populates an even weight 
superposition of the same two single particle states. The 
restriction to only two finally occupied states should be 
justified by a more elaborate treatment of the many-body 
Hamiltonian, and in Sec. Ill, we model the actual process 
in which the atoms are not put directly into two states 
extending over vast and partly non-overlapping regions 
of the atomic single particle quantum states, but where 
the molecular dissociation enforces a localized common 
origin of the matter waves. The problem is solved exactly 
by a Bogoliubov tranformation, and we identify different 
regimes and the validity of the few-mode Ansatz. 



II. TWO-MODE ANALYSIS 

For simplicity we ignore interaction between atoms, 
between molecules, and between atoms and molecules. 
We also make the assumption that only one quantum 
state for the molecules and two quantum states for the 
atoms are relevant. The molecule state should be the 
macroscopically occupied state of the molecular conden- 
sate. The two atomic states are assumed to be se- 
lected by a resonance condition: They could for exam- 
ple be the two lowest 2D single particle eigcnfunctions 
in a harmonic trap [cf. Eq. (Ill below], <&oo(£,y) oc 



exp(-p 2 /2al sc ),<f> i(x,y) cx (F+ iy) exp(-p 2 /2a 2 osc ), 
where p 2 = x 2 + y 2 , and where the oscillator length 
a osc = \J /i/TOatmU; is defined in terms of the atomic mass 
TO atm and the trap frequency w. The lowest state $ o 
has angular momentum m — 0, while <I>oi has m' = 1 
so that they fulfill the condition for angular momentum 
conservation when molecules with M — 1 are dissociated, 
m + mf = M. 

If we can drive the system exactly at resonance, i.e., 
either the atomic and molecular states are degenerate, or 
in a laser induced dissociation process, the field frequen- 
cies involved exactly fulfill the Bohr frequency condition 
for quantum transitions, our model Hamiltonian becomes 



H — (3 ca\h\ + (3* ffaoax. 



(1) 



Here j3 quantifies the strength of the microscopic cou- 
pling, the operator c(&) annihilates (creates) a molecule, 
and the operators So, Si(aJ, h\) annihilate (create) atoms 
in the two atomic states. This Hamiltonian is well known 
in quantum optics where it describes the non-degenerate 
optical parametric oscillator. In general it leads to a com- 
plicated entanglement between the pump beam (c) and 
the signal and idler beams clq, hi which can of course be 
made subject to detailed investigation [5], but in the limit 
of a strong pump, the depletion can often be ignored and 
the pump operators can be replaced by c- numbers. We 
will assume here that the molecular condensate is suffi- 
ciently large that depletion can be neglected, and hence 



we shall consider the simpler Hamiltonian 
H = xa\a\ + x*aoai 



(2) 



The coupling strength \ now includes the c-number de- 
scribing the molecular field. Note that the quadratic 
Hamiltonian ^ does not conserve the number of atoms, 
and the final state appears to be a coherent superposi- 
tion of states with different numbers of atoms. This is, 
however, only because we omit the meticulous writing of 
the associated molecular components of the states, which 
precisely account for the conservation of atom numbers. 
As we shall only be interested in number conserving 
observables such as the atomic density distribution, we 
shall make no errors in applying the symmetry breaking 
Hamiltonian 

Starting in a state with no atoms, the atomic vacuum 
state, in both modes, the time evolution leads to the 
production of the two-mode state 
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where 



-i—rtanhdxIV^)- 

Ixl 



(3) 



(4) 



As discussed in detail in the Introduction, this state 
does not follow from a mean field analysis and, indeed, 
the mean values of the field operators hi vanish exactly, 
whereas the state has a mean atom number of 

(N t ) = (alai) = \s\ 2 /(l - \s\ 2 ) = sinh 2 (\ X \t/h) (5) 

in each mode. The atom number distributions are expo- 
nential (thermal) and thus the fluctuations in the atom 
number are large, Var(iVj) = (Ni) 2 + (N%). 

A number of papers have discussed observational dif- 
ference between systems populating several single parti- 
cle states macroscopically and a coherent superposition 
of the same states. It has been argued [TTJ [T2 |TS] 
that such differences are small or insignificant for systems 
with many atoms. The general argument was supple- 
mented by simulations of actual detection records, where 
the back action due to local measurements on the system 
turned out precisely to establish a definite relative phase 
of the two components. 

The restriction to only two modes implies that the one- 
body density matrix is fully characterized by the 2x2 ma- 
trix with pij — (a\cij). In particular the spatial density 
is given by 



n(x,y) = paa\<&m{x,y)\ 2 + pn\<bm(x,y)Y 
+ 2Re[p m <5> m {x,y)<5>* m {x 1 y)]. 



(6) 



Registration of a particle at position (xd,yd), chosen ac- 
cording to this probability distribution, causes the ap- 
plication of the annihilation operator &oo(xd, 2/d)«o + 
&oi(xd,yd)hi on the state vector expanded in the two 
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mode number components as in (3), followed by a renor- 
malization. We are thus able to compute the up-dated 
2x2 density matrix and by repeating these steps to sim- 
ulate the subsequent detection of a number k of atoms. 
The present problem differs in two way from the simula- 
tions reported in Refs. [TOl [12]: the two modes populated 
are not plane waves but they have different spatial de- 
pendencies, implying that some atomic detection events 
can be ascribed solely to one component (e.g., only the 
m = component contributes to the atomic density on 
the vortex axis) , and the inital state does not have a well 
defined number of atoms, but much larger than Poisso- 
nian fluctuations, and hence the rigid mathematical anal- 
ysis of the emergence of a relative phase in Ref. [H] does 
not apply. We have carried out simulations, and despite 
these two differences, the 2x2 density matrix indeed ap- 
proaches a pure state projector and we obtain detection 
patterns that are compatible with the density of a single 
coherent quantum state. 

III. MULTI-MODE ANALYSIS 

Let us now analyze the problem without making the 
two-mode simplification. Retaining the approximation 
of replacing molecule operators with c- numbers, a full 
many-mode version of Eq.Q reads 

H = d 3 r &(v)h(v)$(v) 

f f 1 (?) 

+ I / d\d\' x(r,r')# + (r)$ + (r') +h.c. L . 

The single particle part consists of kinetic energy and 
external trapping potential 

h(r) = + F cx (r) - A baro . (8) 

The energy offset Abaro is half of the two-atom detuning: 
the energy of a trapped molecule in the vortex state is 
the energy of two free atoms plus 2Abarc- 

Typical traps are well approximated by a harmonic po- 
tential Vcx(r) — rn atm (w^x 2 + uj 2 y 2 + uj 2 z 2 )/2. We will 
consider situations with axial symmetry and we therefore 
let For the purpose of discussing the dis- 

sociation of a molecular condensate with a vortex along 
2, it is useful to make one of two simplifying assumptions 
about the dynamics along z: 

• A quasi cylindrical situation with oj z <C ujj_ such 
that the z component of the momentum, p z , be- 
comes approximately conserved. 

• A quasi two-dimensional situation where lj z ^> uj± 
so that the dynamics along z is effectively frozen to 
take place in a single quantum state. 

We focus on the second case in the following. The modifi- 
cations needed to treat the first case are briefly discussed 
in the Appendix. 



The Hamiltonian ([7| can describe many processes 
where pairs of particles are created: one must choose 
the correct coupling kernel x( r i r ') for the problem un- 
der consideration. When dealing with dissociation of 
molecules, ,\(r, r ') depends on the wavefunction of in- 
dividual molecules (both relative motion and center-of- 
mass) , and on the possible spatial variation of the exter- 
nal fields mediating the dissociation. We will assume the 
relative motion molecular wavefunction to be very well 
localized compared to other length scales in the problem. 
This means that the dependence of x(r, r') on the separa- 
tion r — r' can be approximated by a delta-function. We 
further assume that the center-of-mass molecular wave- 
function and the relevant external fields are rotationally 
symmetric around the z-axis so that 

X (r,r')-^(r-r')Xc m (p,z)e lM *, (9) 

in usual cylindrical coordinates where x = p cos <f) and 
y = psin0. In the simplest case the dissociating external 
fields are spatially uniform over the extend of the molec- 
ular condensate. Then M is the charge of the molecu- 
lar vortex state and the remaining spatial variation of 
Xcm(p, z) is simply given by the norm of the c- number 
field that describes the molecular condensate (squareroot 
of the density) [20] . To specify the overall strength of the 
coupling a microscopic model of the dissociation must be 
decided upon. In Ref. [2], a two-photon Raman process 
between a bound molecular state and "free" (except for 
external trapping) atoms is considered. The dissociation 
is via an excited molecular state, which should ideally be 
only negligibly populated since spontaneous decay from 
it will be an unwanted loss mechanism. The particu- 
lar example in Ref. |14j shows that a peak value of our 
X comparable to the trap level spacing is fully realistic. 
This is the regime we focus on below. 

We emphasize that this physical modelling does not fa- 
vor any particular pair of atomic modes, and the Hamil- 
tonian ^ cannot be brought on the form ^ . Atoms are 
coherently prepared from the entire region populated by 
the molecules, but pairs are initially much more tightly 
located both radially and azimuthally than the single 
mode functions suggested in the previous section. In 
particular, the atoms are not restricted to low values of 
the angular momentum quantum number. The tempo- 
ral evolution of the system due to the kinetic energy and 
external potential, and the role of energy conservation, 
however, favors the population of only few modes as we 
shall see below. 



A. Frozen z dynamics 

In the case of tight z confinement for the atoms, we 
get a simple description. Let the only accessible atomic 
z mode be 4> z , J dz\(f> z \ 2 = 1. We can then expand the 
atomic field operators on the discrete set of mode opera- 
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tors 

*(r) = £ a nm x $ nTO (p) x J^e™* x <j>,(z) (10) 

nm 

where the radial modefunction Q nrn are 

/ ri I 



2(n + |m|)!et. 



2 

osc 



x exp 



IP 2 



„ \ M / „2 

— ) 4r' (4- 



(11) 



with L™ the m'te associated Laguerre polynomum of or- 
der n and a osc — -\/fi/m a t m w is the oscillator length. The 
commutation relations of the a nm are 



= 



(12) 



The Hamiltonian becomes 



}, (13) 



where K nn i m is defined via 



A' 



(kdpp Xcm(p, z) \<l>(z)\ 2 $* m (p)$ n ' M -m(p)- 

and £7 nm is given by 

£„ m = (2n + |m| + 1) Hu - A. 



(14) 



(15) 



Here A is the effective detuning, adjusted for the energy 
associated with <b z : 



A = Abarc - / dz 



P 



1 



2m + 2^ ) +" (16) 



IV. BOGOLIUBOV DIAGONALIZATION 



The Hamiltonian ( 13 1 is a quadratic form of creation 



and annihilation operators and as such it can in gen- 
eral be decoupled to a collection of independent harmonic 
oscillators by a Bogoliubov transformation. Please note 
that this transformation is usually applied in connection 
with the Bogoliubov approximation in theoretical stud- 
ies of Bosc- Einstein condensed systems. The Bogoliubov 
approximation provides a quadratic Hamiltonian by ap- 
peal to the macroscopic population of the condensate 
mode. Here the quadratic form of the Hamiltonian is 
rather a consequence of (i) that two atoms are created 
in each fundamental dissociation process and (ii) that we 



ignore the dynamics (in particular the depletion) of the 
molecules and describe them by a stationary c-number 
field. One could in principle extend our approach to also 
treat the interaction among the created atoms in a Bo- 
goliubov approximation, but since we focus on situations 
with initially no atoms present this would neither be nec- 
essary nor indeed well justified. For an example of time- 
dependent quadratic Hamiltonians, sec the work by Zih 
et al. on colliding condensates [TBI 116) . 



A. Two-mode Bogoliubov 

Let us first look at the simple case of two modes like in 
Sec. [TTJ In addition to the pair creation and annihilation 
terms of Eq. ^ we include mode energies Hu>o and hu>\: 



H = hojQ ( CIqCIq 



+ X%a\ + x*aoai. 



huj\ I a\a\ 



(17) 



A useful and compact notation is 



1 



H = -A^hA 
2 



(18) 



with 



A = 



At = [ 



a\ a\ a oi] , 



(19) 



and 



Hluq x 

tkJi x 

x* 7^0 

.x* o o hu>i. 



(20) 



We now seek to simplify H by a Bogoliubov transforma- 
tion, i.e. we define new creation and annihilation opera- 
tors implicitly by 



U V* 
V U* 



(21) 



where the 2x2 matrices U and V should fulfill 

U*U-V*V = h 
U T V - V T U = 



(22) 



in order for bo and b\ to have the commutation rela- 
tions for independent bosonic creation and annihilation 
operators [see Eq. ( 12 1] . Note that in general the Bogoli- 
ubov transform is not simply a choice of different spatial 
modes. 
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1. Detuning dominated case 2\\\ < \hcoo + hu>i\ 

Depending on the strength of the coupling \\\ relative 
to the two-atom detuning \hcoo + fiuii\ the Hamiltonian 
can be written in one of two standard forms. If 2\\\ < 
\hLUQ+huji | , H can be written as a sum of two independent 
oscillators 



H = Ac I blk 



(23) 



with 



A ,i = ± -(hojo ~ f^Ji) 



4|xl : 



(24) 



I Huq + huii I 



The Bogoliubov transformation is given by: 

U = 



cosh r 
cosh r 



and 



V 



T^y sinhr 



=F fr sinh r 
T Ixl 







(25) 



(26) 



where the upper (lower) sign should be chosen for Huq + 
hu>i positive (negative). The squeezing parameter r is 
defined through 



tanh r 



2|X| 



\hiv + fax | + ^\hui + hwi\ 2 - 4\x\ 2 



(27) 



We see that as 2\x\ approaches \hu>o + Huii\, the coeffi- 
cients in the transformation diverge, i.e. the decoupled 
modes become infinitely squeezed. 

Because of the diagonalized form of the Hamiltonian 
( 23 1 , the time evolution will be independent for the two 
Bogoliubov modes: Each will simply behave as a har- 
monic oscillator. This means that the dynamics of all 
quantities will be oscillatory with two fundamental fre- 
quencies given by Ao,i/?L 

Note, that even if both mode energies Hu>o,i are posi- 
tive so that a pair of atoms actually has a higher energy 
than a molecule, one of the eigenvalues Xi can become 
negative. This happens if \x\ 2 > hio^huii and in that 
case the system will be thermodynamically unstable as 
it is energetically favourable to create the corresponding 
kind of quasi-particles. 



real positive coefficient) into two symmetrically detuned 
modes: 



H = Re 
with 



[A] (blh - blk) + Im[A] (P b\ + boh) (28) 



A =-(hu>x - huj ) 



(29) 



The transformation is still of the form given in Eqs. (25 1 
and ( |26[ ), but now the squeezing parameter r should be 
found from: 



tanh r 



\hujQ + hu>i\ 



2\ X \ + v/4|x| 2 - l^o + ^i I 



(30) 



We see that as before, r diverges when 2|x| approaches 
\hujQ + hwi\. At exact two-atom resonance, i.e. Hloq + 
hu>i = 0, the Bogoliubov operators are simply the original 
So and &i. 



The transformed Hamiltonian ( 28 1 consist of two com 



muting terms and the first one even conserves the num- 
ber of quasi-particles Nf, = b^bo + b\b\. The second term 
leads to unbound creation of quasi-particle pairs, in fact 
(Nb) will grow as sinh 2 (Im[A]<). This production of quasi- 
particles translates to a similarly unbounded growth in 
the number of atoms (until depletion of the molecular 
condensate renders our simple Hamiltonian invalid) , and 
because of the Bogoliubov tranformation, the atom num- 
bers will show oscillations around their general growth. 



B. Multi-mode Bogoliubov 

The generalization to many modes is relatively 
straightforward. By Bogoliubov transformations of the 
original set of creation and annihilation operators, H can 
be written as a sum of a number of independent oscil- 
lators and a number of pairs of symmetrically detuned 
modes into which quasi-particles are created: 

+ Re t A f] (Hah* - (31) 

p 

+ Im[A p ] (P pl P p2 + b pl b p2 ) . 



2. Coupling dominated case 2\x\ > \hcuo + hu>i\ 

If 2 Ixl > \tuJo + the simplest form attainable by 
Bogoliubov transformations involves pair-creation (with 



The dynamics in the decoupled modes is oscillatory while 
ever more quasi-particles will be created in the paired 
modes. Note that it is perfectly possible that all modes 
can be decoupled so that the unitary dynamics is purely 
oscillatory. 
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C. Consequences of the azimuthal symmetry 

We consider a situation with axial symmetry: The sin- 
gle particle part of the Hamiltonian, h(r) of Eq. (J8j, is 
invariant under rotations around the z-axis. This leads 
to the usual block-diagonal form of the single particle 
Hamiltonian with respect to the quantum number to, 
the eigenvalue of the ^-component of the angular mo- 
mentum operator. At the same time, the molecular field 
[x of Eq. ([9])] changes simply by a phase factor exp(iM(f>) 
when the system is rotated by an angle <f>. This leads to 
an exact selection rule on the pairs of modes that are pop- 
ulated in the dissociation process, namely m + to' = M. 
For odd M this imposes an "off-diagonal" form on the 
pair creation, as the atoms created by a dissociation pro- 
cess must necessarily end up with different integer quan- 
tum numbers, to and M — to, whereas for even M the 
two atoms can both end up in an m — M jl mode. 

In summary, each to subspace is coupled only to itself 
by the single particle part of the Hamiltonian and the 
dissociation can either couple it exclusively to itself or 
exclusively to one other to subspace. In this way the full 
problem splits into a series of smaller problems consisting 
of either one or two TO-subspaces. 

Now consider the problem of two coupled subspaces. 
As pairs are always created with an atom in each sub- 
space, there is an obvious conservation of the difference 
in the total number of atoms in the two subspaces. This 
can also be phrased as the invariance of the evolution 
under shifts in the relative energy of the two subspaces 
as long as their total energy is maintained. Among other 
consequences, this implies that there are no first-order 
coherences between the two subspaces and that the one- 
body density operator is correspondingly block-diagonal. 
It is, in fact, possible to show an even stronger, dynami- 
cal symmetry in the evolution: If the two subspaces are 
initially unpopulated, they will at all times remain uni- 
tarily equivalent. In particular, their one-body density 
matrices will have identical spectra of eigenvalues 17J. 



V. RESULTS 

As discussed above, once the Hamiltonian has been 
brought to a standard form by a Bogoliubov transforma- 
tion, a number of properties of the system are immedi- 
ately clear. In particular, the oscillatory or unbounded 
behaviour depends on whether all modes can be decou- 
pled. In Fig. [I] we show the real and imaginary parts 
of the Aj's as a function of the detuning A. When A 
is increased, all atomic modes move linearly down in en- 
ergy according to the single atom Hamiltonian (7). As a 
pair of modes with m + m' = M = 1 gets close to the 
two-atom resonance, i.e. to having equal and opposite 
energies, the coupling has the effect of bending the quasi- 
particlc energy levels further downwards. The modes be- 
come increasingly squeezed but remain uncoupled - they 
are in the detuning dominated regime of Sec. IV A 1| At 




FIG. 1: Bogoliubov eigenvalues as functions of the detuning 
A for the Hamiltonian ( 13 1 with a coupling field Q with 
M = 1. For concreteness, we assume a strength and a ra- 
dial dependence of the coupling such that j ' dzXcm\4>z\ 2 = 
2huipexp(—p 2 /a,a ac )/a OB c- The lines show Re(A), which are 
the mode energies while, in the coupling dominated regime, 
the shaded areas around them have a height of Im(A) to in- 
dicate the coupling between pairs of modes. The darkness 
of the shading indicates the m-values in the pairs of cou- 
pled modes: the darkest shading is for m = 0, m' = 1 pairs, 
m = — l,m' = 2 pairs are a tone brighter, etc. Note that 
members of coupled pairs he symmetrically around A = 0. 



the point of perfect two-quasi-particle resonance, it is no 
longer possible to decouple the two modes and the sys- 
tem moves into a regime with quasi-particle pair produc- 
tion, the coupling dominated case of Sec. |IV A 2| When 
a pair of quasi-particle modes become unstable, it is sig- 
nalled by the corresponding Aj's attaining a finite imag- 
inary part. In Fig. [I] this is plotted as a shaded region 
around the corresponding energy curves. The height of 
the shaded region signifies the magnitude of Im[A]. As A 
is further increased, the pair production strength also in- 
creases. It has its maximum approximately at the point 
of two-aiom resonance, i.e., when the two bare atomic 
modes becomes resonant with the molecule energy. 

As an example, let us follow the lowermost curve 
through the figure. At the left, at low values of A, 
the curve corresponds to the atomic harmonic oscilla- 
tor ground state in the trap with n = and m = 0. 
When A gets close to 1, pair production together with 
the (n = 0,m = 1) state (second lowest curve) becomes 
resonant. Note that the (n = 0,m — 1) state is degener- 
ate with the (n = 0, to = — 1) for very low A. At approx- 
imately A = l.bfuu [i.e. £00 + E i = 0, cf. Eq. ( fl5| ] the 
pair production has its maximum. A small shift from the 
simple estimate is due to the multi-mode character of the 
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problem. The next resonance for the (n = 0, m = 0) level 
is centered around at A = l.hhw [i.e. E 00 + En = 0]. 
Note that around A = 2.5/kj, the production of pairs 
with (n — 0, m = —1) and (n = l,m = 2) is also res- 
onant. The third resonance, where Eqq + E21 = 0, is 
at A = 3.5/kj. Here several other resonances are also 
present. 

In Fig. [2] we plot the total number of atoms and in 
Fig.[3]the fractional distribution on single-particle modes 
as a function of time for three different values of A. For 
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FIG. 2: (color online) Total number of atoms as a function of 
time for three values of A. The coupling is as in Fig. [I] 

A = 0.80/kj the dynamics is detuning dominated and os- 
cillatory. All modes can be decoupled and the resulting 
quasi-particle modes are only slightly squeezed resulting 
in a very modest production of atoms. This almost van- 
ishing population mainly occupies two modes: one in the 
m = manifold and one in the m' = 1 manifold. These 
modes constitute the pair that is closest to resonance 
and their oscillation period is of the same order of mag- 
nitude as the trap period. The fractional occupation of 
these modes is close to 50% except for times close to 
the minima in their oscillation (e.g. Lot ~ 7 here). At 
A = lAOtuo, the system can no longer be decoupled: A 
single m — 0, m! — 1 pair of modes has just become 
unstable and real pair production takes place. Therefore 
this pair of modes quickly becomes dominant. Finally, 
for A = 2.\bh\o the evolution is more complex. From 
Fig.[l]one can tell that there will be real pair production 
into a pair of modes with m = — 1, m! = 2. However, a 
pair of modes with m — 0, m! — 1 also becomes unstable 
if A is increased just a little bit. It turns out that the 
m = 0,m' = 1 pair dominates at early times, but later 
most atoms are produced in the m = — 1, m! = 2 pair of 
states. 

The above results all concern one-body properties. An 




FIG. 3: (color online) Condensate fraction, i.e. relative weigth 
of the largest eigenvalue of the one-body density operator 
compared to the sum of all eigenvalues (total number of par- 
ticles. Results are plotted for three different values of A. As 
discussed in Sec. |IV C| the conservation of angular momen- 
tum dictates a block-diagonal form of the one-body density 
operator and that pairs of blocks for which m + m' = M will 
have identical spectra. Fully drawn lines correspond to m = 1 
and m! = 0, dashed lines to m = 2 and m' = —1, and dotted 
lines to m = 3 and m' — —2. 



illustrative way to quantify some of the two-body cor- 
relations in the system is to calculate the conditional 
density distribution, i.e., the density distribution result- 
ing after the detection of a single atom at some position 
( x d,yd)- In Fig. [1] we plot the result of such a calcu- 
lation for A = 2.15huj and tut = 2,4,6. The upper 
row shows the normalized density distribution, that is, 
the probability distribution for the detection of the first 
atom. In the lower row, an atom has been detected at 
(xd,Ud) — (l,0)a osc and the resulting conditional prob- 
ability distribution for the next detection is shown. For 
all three times, this conditional distribution is dramat- 
ically different from the original one. Close to (xd,yd), 
there is generally a significantly higher probabilty den- 
sity, reminiscent of the familiar (thermal state) bunching 
of bosons. However, when higher angular momentum 
states are populated, additional peaks appear as the de- 
tection induces coherences between more m values. From 
a calculation of the full one-body density operator we can 
also find the new "condensate fraction" and for all three 




FIG. 4: (color online) Detection probability distribution (normalized density distribution) before (upper row) and after (lower 
row) detection of a single atom in (xd,Vd) = (l,0)a OS c- The detuning A was chosen to be 2.15hu>. The three columns show 
results where the system was allowed to evolve for a time u>t = 2, 4, and 8, respectively. Other parameters were like in the 
previous figures. Before the detection of the first atom, the density distribution is rotationally symmetric and it is equally likely 
to make the first detection at any azimuthal angle, <f>. According to Fig. [3] the one-body density matrix at ut = 2 is an almost 
equal mixture of m — and m! = 1 atoms. In each of these m-manifolds, a single radial mode dominates and accounts for 
35% of the total number of atoms. In the lower panel, we assume that an atom has been detected with cj> — 0, more precisely 
at (xd,yd) = (l,0)a osc . The conditional probability distribution for the next detection has become asymmetric indicating a 
significantly increased probability to detect the next atom close to the first one. Diagonalization shows that the one-body 
density matrix now has a "condensate fraction" of almost 60%. At cut = 4 a significant number of m = — l,m' = 2 pairs 
have been created. The density before detection has a pronounced dip at the center due to centrifugal effects. The conditional 
density after the detection is rather complex: it is in fact not even symmetric around the t/-axis. Because of the larger number 
of modes involved, the condensate fraction is only 43%. Finally, at cut = 8 the m — —l,m' = 2 pairs dominate. After the 
detection, three almost equally strong peaks are seen, signaling the appearence of coherences between two modes with angular 
momentum differing by 3 units. The condensate fraction after the detection is now again almost 60%. 



cases this number is significantly increased compared to 
the situation before the detection (cf . caption of Fig. [4]) . 
This is similar to what has been proposed in the case of 
number state condensates [T5] , 



VI. CONCLUSION 

We have presented a qualitative discussion and a quan- 
titative analysis of the break-up of a molecular con- 
densate into an atomic system, in the case where the 
molecules occupy a vortex state with unit angular mo- 
mentum. The rotational symmetry of the problem causes 
the fragmentation of the atomic system into (at least) two 
spatial components, and our multi-mode analysis quan- 
tified the dynamics of these components. Particle detec- 
tion at random locations will generally break the rota- 
tional symmetry of the state, and both our simple two- 
mode model and our more general analysis show that 
such detections will indeed cause the system to approach 
a single macroscopically populated quantum state with a 
preferred angular dependence. 



Several interesting possibilities exist for further stud- 
ies. A more technical issue concerns our assumption of 
non-interacting atoms. As shown in [7] interactions have 
non-trivial consequences for the quantum correlations in 
the system. Depletion of the molecular condensate [14] 
affects the quantum statistics of the atomic state, and 
moreover, the explicit inclusion of the molecular field 
will also enable processes where atoms in different frag- 
ments recombine and form molecules with new angular 
momenta, respecting always the overall rotational sym- 
metry of the system, but introducing multi-atom correla- 
tions in the system. Within our simple model with non- 
interacting atoms and undepleted molecular systems, it 
is a natural question to ask, what atomic states will result 
from the dissociation of a molecular vortex lattice. Will 
two fragments be sufficient to describe such states? We 
assumed a cylindrical symmetry and referred extensively 
to angular momentum conservation in the arguments of 
this paper, but in an interacting system vortices also ex- 
ist under non-symmetric confinement, and our symmetry 
argument should be adequately modified into an argu- 
ment referring to the phase topology. Presumably one 
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would still see the formation of a fragmented system with 
components with and without vorticity, but the need for 
interactions among the atoms for the stability of such a 
state is an interesting issue. 
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Under the above assumptions, the Hamiltonian ^ then 
simplifies to the form: 



t'nmk 



nmk 



(M-m)(-fc) 



mk nn' 



+ h.c.}, (A.3) 



APPENDIX: CYLINDRICAL SYMMETRY 



where 



The case of strict cylindrical symmetry can be treated 
in a very similar fashion as the one of frozen z dynamics 
of Sec. |III A| Assume that potential is axially symmetric 
and flat along z: oj z = 0. The coupling is also assumed 
to be z independent, \ cm (p,z) — K{p). We then use a 
quantization box of length 2L along z for the atoms. The 
field operator is expanded as follows 



oo oo 



n— m— — oo k— — oo 



} i7rkz/ L 



(A.l) 



where the radial modefunctions where given in Eq. (Ill 



The mode creation and annihilation operators satisfy 
usual commutation relations 



Qnmk-) Q"n'm'k' 



5 n'm'k' 



= o 

— $ fin' $mm' $kk' ■ 



(A.2) 



E nmk = 2n + \m\ + 1 + ~ (~) k 2 - A 



(A.4) 



and 



A" 



dzdppK(p) — $* nm (p)$ n , M _ m (p). 



(A.5) 



The main difference to the case of a single active z 
mode treated above is that excess molecular energy can 
now be transformed to z kinetic energy of the atoms. 
Thus the picture of isolated resonances breaks down, even 
at low coupling strengths. Many atomic modes will par- 
ticipate in the dynamics and it will take more detections 
to build up a sizable condensate fraction. 
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